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Summary 

Results are presented of elastic, elastic-plastic, and elastic- 
plastic-creep analyses of a test-rig component of an actively 
cooled cowl lip. A cowl lip is part of the leading edge of an 
engine inlet of proposed hypersonic aircraft and is subject to 
severe thermal loadings and gradients during flight. Values 
of stresses calculated by elastic analysis are well above the 
yield strength of the cowl lip material. Such values are highly 
unrealistic, and thus elastic stress analyses are inappropriate. 
The inelastic (elastic-plastic and elastic-plastic-creep) analyses 
produce more reasonable and acceptable stress and strain 
distributions in the component. Finally, using the results from 
these analyses, predictions are made for the cyclic crack 
initiation life of a cowl lip. A comparison of predicted cyclic 
lives shows the cyclic life prediction from the elastic-plastic- 
creep analysis to be the lowest and, hence, most realistic. 

Introduction 

A significant research effort was launched recently at NASA 
Lewis Research Center to develop state-of-the-art technologies 
for hypersonic flight between Mach 3 and 25. At such high 
speeds, the leading edges of hypersonic aircraft are subjected 
to high heat fluxes and, thus, high temperatures and severe 
thermal gradients. To achieve high inlet aerodynamic 
performance, not only must the high heating rates be tolerated, 
but also distortions caused by thermal warping of the structure 
must be minimized. Consequently, the need arises for the 


development of actively cooled leading edges fabricated from 
specialized materials with innovative cooling concepts to 
enable the structure to withstand the severe environmental 
conditions. The details of different cooling concepts proposed 
and studied under a NASA Lewis sponsored program, called 
COLT (Cowl Lip Technology Program), were presented by 
Melis and Gladden. 1 

Under the severe thermal loadings and gradients experienced 
by the component, a propensity for significant inelastic 
deformation of the material exists. The intent of the present 
work was to perform linear (elastic) and nonlinear (elastic- 
plastic and elastic-plastic-creep) structural analyses of the cowl 
lip for a crossflow cooling concept, 1 and to assess potential 
durability to low-cycle fatigue cracking. In crossflow cooling, 
the direction of coolant flow is perpendicular to that of the 
hot gas. (See Figure 1.) A comparison of service lives 
predicted from the outputs of the three structural analyses is 
also presented. The three-dimensional finite element model 
and steady-state heat transfer results were obtained from 
Reference 1 . The imposed thermal loadings obtained from a 
rig test of a cowl lip are shown in Figure 2. Measured gas 
and wall temperatures are shown for one complete cycle. For 
ease of calculations, the leading-edge metal temperature was 
approximated by a simplified temperature cycle. This sim- 
ulated thermal loading cycle, shown in Figure 2(c), was used 
in subsequent calculations. Nonlinear variations in the material 
constants with temperature were accounted for in the 
computations. The life assessment calculations for the cowl 
lip were performed by using data from the existing 
literature. 2 Details of the finite element model and thermal, 
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structural, and life analyses of the cowl lip are presented in 
the following sections. 

Analyses 

Finite Element Model 

A three-dimensional, finite element model of the cowl lip 
(Figure 1) was constructed by Melis and Gladden. 1 The 
model consists of 3294 eight-noded, solid, isoparametric 
elements and 4760 nodes. The dimensions of the rig component 
are 6 by 1.5 by 0.25 in. (15.2 by 3.8 by 0.6 cm). The cowl 
lip finite element model consists of the central 2-in. (5-cm) 
portion of the component. This small section was analyzed 
to avoid the difficult-to-quantify constraining effects of the 
ends. A considerably large number of elements are required 
to handle the severe thermal loadings imposed on the com- 
ponent. The severity of thermal loading, together with the fine- 
grid finite element model, makes the nonlinear structural 
analyses computationally intensive. However, taking advantage 
of the symmetry along the thickness reduces the computation 
time appreciably. 

Heat Transfer Analysis 

Reference 1 describes a boundary layer analysis conducted 
to compute the film coefficients. The authors used the 
boundary-layer heat transfer code STAN5 3 to estimate the 
film coefficients on the hot-gas side of the cowl lip. Correlative 
techniques 4 ’ 5 were employed to determine coefficients on the 
coolant side. The stagnation points and leading-edge-region 
film coefficients were ascertained from a “cylinder in 
crossflow” correlation. Using the film coefficients determined 
by these procedures, the authors performed a steady-state heat 
transfer analysis for the cowl lip. An instrumented rig 
component was tested in the Hot Gas Test Facility at NASA 
Lewis Research Center. 6 Thermal loadings are displayed in 
Figures 2(a) and (b). 

Structural Analyses 

Linear elastic and nonlinear elastic-plastic analyses were 
performed by using the finite element program MARC. 7 The 
material of the cowl lip was copper. The temperature 
dependence of Young’s modulus E and the coefficient of 
instantaneous thermal expansion a, taken from Freed and 
Verrilli, 8 are listed in table I. The elastic Poisson’s ratio v 
was assumed to be 0.34 for all temperatures. 

The structural response of the cowl lip was calculated by 
using an eight-noded solid brick element. This element was 
compatible with the element used for the heat transfer analysis. 
The steady-state temperature values (between 0.75 and 
2.25 sec in fig. 2(c)) were obtained from Reference 1. A linear 
interpolation technique was then employed to obtain the 
transient temperature values (in the periods between 0 and 


TABLE I.— MATERIAL CONSTANTS FOR COPPER 


Young’s modulus, E , ksi a 

18 504.0 - 3. 164 P> - 0.0022 r* 

Coefficient of thermal expansion, 

a 1.59x10-5 + 2.78X10- 9 T 

Poisson’s ratio, v 

0 34 

A, sec~' 

3 0X10 24 

n 

5 o 

m 

0 35 

Activation energy, J/mole 


Lattice diffusion, Q L 

2.0x105 

Pipe diffusion, Q ? 

1.2x105 

X 

1 8xl0 -5 

Boltzmann’s constant, J/mole K .. 



a 1 ksi = 6.895 MPa. 

b T is temperature in degrees Celsius. 


0.75 sec and 2.25 and 3 sec) to describe the temperature 
distribution over the complete thermal loading cycle. For the 
nonlinear elastic-plastic analysis, a combined isotropic and 
kinematic hardening plasticity law was employed. The tem- 
perature and hardening effects were incorporated in the analysis 
through the input for the finite element program MARC. 

To investigate the time-dependent response of the cowl lip, 
a nonlinear combined elastic-plastic-creep stress analysis was 
performed. The following creep law was used in the analysis: 

e = a) 

where 

" exp (jf) + * exp (if) <2) 

In these equations, o e is the effective stress; E is Young’s 
modulus; t is the time; A, n and m are material constants; Q L 
and Q P are the activation energies for lattice and pipe 
diffusion, respectively; X is the constant to account for the 
activation energy for pipe diffusion that must be included in 
the temperature function 6 above a homologous temperature 
of about 0.5 (see Figure 3); k is Boltzmann’s constant; and 
T k is the temperature in degrees Kelvin. This form of the 
creep law was obtained by using the experimental data of 
Barrett and Sherby 9 and Freed and Verrilli 8 and is shown in 
figures 3 and 4, respectively. The material constants appearing 
in the creep law for copper are listed in table I. 

Elastic analysis.— Figures 5 and 6 depict, respectively, the 
elastically calculated stress and total strain distributions along 
a cross section of the cowl lip at thermal steady state. The 
total strain in this work does not include the thermal component 
of strain. Only a cross section of the component is shown in 
these figures. Since the largest compressive stresses and total 
(mechanical) strains occur along the leading edge (z-direction), 
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TABLE II. -COMPARISON OF HYSTERESIS LOOP CHARACTERISTICS 


Analysis 

Effective stress 
at maximum 
temperature 

Effective stess 
range 

Total 

mechan- 

ical 

strain 

range 

Elastic 

strain 

Plastic 

strain 

Creep 

strain 

Cycles 

to 

failure, 


ksi 

MPa 

range, 

Ae c , 

range, 

Ae PP 

range, 


ksi 

MPa 

Ae PC 

N r 

Elastic 

-80.4 

-554.4 

105.4 

726.7 

0.0062 

0.0062 



2300 

Elastic- 

plastic 

-1.4 

-9.7 

8.8 a 

60.7 

.0085 

.0005 

0.0080 


1000 

Elastic- 

plastic- 

creep 

-1.3 

-9.0 

8.7 

60.0 

.0092 

.0005 

( — ) .0080 

(=5 ).0007 

800 


a Stress range between points A and B of hysteresis loop of figure 12. 


graphs showing these distributions are included in the dis- 
cussion. The magnitude of the largest stress predicted from 
the elastic analysis (see Figure 5) is much higher than the yield 
strength values for the copper material. This clearly indicates 
the inappropriateness of an elastic stress analysis for this severe 
thermally driven problem. The subsequently large imposed 
mechanical strains at these high temperatures are well above the 
yield strain and hence are expected to result in significant inelastic 
flow, necessitating a nonlinear inelastic stress-strain analysis. 

Elastic-plastic analysis.— The stress and strain results of an 
elastic-plastic analysis of the cowl lip problem are shown in 
Figures 7 and 8. These figures exhibit the stress and strain 
distributions along the leading edge (z-direction). The analysis 
shows the reduction of maximum compressive stress to a more 
reasonable value of 28.5 ksi (196.5 MPa) from about 80.4 
ksi (554.4 MPa) obtained from the elastic analysis. From 
Figure 8, the largest compressive strain occurring at the 
leading edge is greater than that obtained from the elastic 
analysis. This indicates the existence of compressive plastic 
strains along the edge. 

Elastic-plastic-creep analysis .— To study the effects of time 
on stress and strain distributions in the cowl lip, a time- 
dependent creep analysis was also performed. The results are 
displayed in Figures 9 and 10. These figures show that the 
maximum compressive stresses and total strains still occur at 
the leading edge. The magnitude of maximum stress decreases, 
whereas the magnitude of total strain increases. This increase 
in strain magnitude implies the existence of compressive creep 
strains along the leading edge. During the steady state of 
thermal loading (i.e., between 0.75 and 2.25 sec), no further 
redistribution of stress due to creep is observed. This is made 
clear by comparing figures 9 and 1 1 , which show the stress 
distributions in the cowl lip at 0.75 and 2.25 sec, respectively. 


effective* stress and strain from the elastic-plastic-creep 
analyses are calculated for the “critical” cowl lip location (i.e., 
where the total strain range is maximum). In establishing these 
loops, the algebraic signs of effective stress and strain are 
dictated by the signs of the maximum principal stress and 
strain. This is in accordance with the recommendations of 
Manson and Halford 10 for dealing with multiaxial stress- 
strain states. Critical-location stress-strain hysteresis loops for 
the first two complete loading cycles are presented in Figure 12. 
Table II shows a comparison of hysteresis loop characteristics 
(stress range, total strain range, elastic strain range, inelastic 
strain range, plastic strain range, and the approximate creep 
strain range) for the three different (elastic, elastic-plastic, and 
elastic-plastic-creep) structural analyses conducted. The elastic 
analysis calculates the smallest total strain range, whereas the 
elastic-plastic-creep analysis produces the largest. Because the 
latter strain range is about 50 percent larger than the smallest, 
the cowl-lip structural geometry apparently offers only limited 
constraint to the inelastically deformed material at the leading 
edge. Hence, the concept of total strain invariance suggested 
by Manson 11 for highly constrained, thermally dominated 
loading problems should not be used for the limited geometric 
and thermal constraint found in the current problem. 

The difference in inelastic strain range between the elastic- 
plastic and the elastic-plastic-creep results can be attributed 
to creep deformation; the creep strain range listed in the table 
is defined simply as this difference. The large inelastic strain 
range is dominated in this instance by time-independent 
plasticity. Furthermore, because of the very large strain range, 
the loops of figure 12 exhibit a high degree of symmetry about 
the zero stress axis. This gives rise to remarkably well- 
balanced hysteresis loops that exhibit little of the unbalanced 


Life Analysis 

The stress-strain results from the structural analyses were 
used in estimating the cyclic crack initiation life of the cowl 
lip. The thermomechanical fatigue (TMF) hysteresis loops of 


Effective stress a e = — [(ffj - a 2 ) 2 


effective strain e e 



+ (a 2 - o 3 ) 2 + (a 3 - a,) 2 ] 1/2 ; 
( f 2 “ e 3>2 + ( e 3 ~ e i) 2 ] * 
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characteristics (CP- and PC-type f hysteresis loops) normally 
associated with thermomechanical cycling at much smaller 
strain ranges. 12 Consequently, the hysteresis loops are 
composed of inelastic strain ranges of the balanced PP- and 
CC-types. The CC portion of the inelastic strain range is 
imposed quite naturally near the very highest temperature, 
whereas the PP strain range is imposed over the entire 
spectrum of temperatures from lowest to highest. 

Because of the large inelastic strains imposed at the critical 
location, a strain-based, high-temperature, life prediction 
approach was sought. Ideally, the most accurate approach 
would be to generate TMF data for the alloy in question. 
Unfortunately, TMF data are virtually impossible to generate 
in the laboratory under the high loading and heating rates 
encountered in the current problem. Thus, approximations 
have been necessary to establish a relation between strain range 
and cyclic crack initiation life that would be applicable for use 
in predicting life for the current problem. 

Isothermal, strain-controlled, low-cycle fatigue results have 
been reported by Conway, Stentz, and Berling 2 for annealed 
oxygen-free, high-conductivity (OFHC) copper at 1000 °F 
(538 °C) for a nonoxidizing, argon atmosphere. These results 
are shown in Figure 13. At the strain range of current interest 
(approximately 0.009), the frequency of loading is about 
0.1 Hz. This frequency is slightly lower and thus more 
degrading than that of the cowl-lip TMF cycle (0.33 Hz, 
nonoxidizing, gaseous hydrogen environment). Because of the 
relative temperature insensitivity of the tensile ductility of 
OFHC copper and dilute alloys of copper such as zirconium- 
copper (variations measured at NASA Lewis range from 81 
to 89 percent reduction of area over the temperature range 
of concern), low-cycle fatigue resistance is also expected 14 
to be insensitive to temperature over the range of current 
interest. Thus, the fatigue curve of Figure 13 serves as a 
reasonable approximation to the TMF resistance of the cowl 
lip problem at hand. Respective total strain ranges of 0.0062, 
0.0085, and 0.0092 for the elastic, elastic-plastic, and elastic- 
plastic-creep analyses correspond to cyclic crack initiation lives 
of 2300, 1000, and 800. Since the elastic-plastic-creep analysis 
is judged to give the most realistic structural analysis results, 
the corresponding life of 800 cycles to failure is also judged 
to be the most realistic estimate of expected lifetime. Note that, 
for the present problem, use of an elastic analysis to obtain 
the total strain range (i.e., using the concept of total strain 
invariance 11 ) would have resulted in a life prediction of about 
a factor of 3 too high (i.e., unconservative). It should be noted 
that the predicted life of 800 cycles is based solely on a low- 


iThe CP, PC, CC, and PP terminology for describing the hysteresis loops 
derives from the strainrange partitioning life prediction method of Manson, 
Halford, and Hirschberg. 13 The letter P refers to time-independent plasticity 
and C to time-dependent, thermally activated, creep deformation; the first 
letter refers to the tensile and the second to the compressive inelastic strains. 


cycle fatigue criterion for failure. Other failure modes such 
as thermal ratcheting could further degrade cyclic durability. 

Conclusions 

Elastic, elastic-plastic, and elastic-plastic-creep analyses for 
the cowl lip were performed. A purely elastic stress-strain 
analysis for the problem is found inappropriate, since it 
calculates values of maximum stresses along the leading edge 
that are much larger than the yield strength of the copper used 
to manufacture the cowl lip. The elastic-plastic or elastic- 
plastic-creep analyses are found to produce more reasonable 
and hence acceptable values of stresses. This clearly indicates 
the necessity for performing inelastic analyses for severely 
loaded structural components. 

The estimations for cyclic crack initiation life of the cowl 
lip were made by using the results from the elastic, elastic- 
plastic, and elastic-plastic-creep analyses. The lowest expected 
life of 800 cycles is obtained by using the results from the 
most realistic elastic-plastic-creep analysis. Use of an elastic 
analysis gives an expected life of 2300 cycles, about a factor 
of 3 higher than that obtained from the elastic-plastic-creep 
analysis. 
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Figure 1.— Cowl lip Finite element model. 



(a) MEASURED HOT GAS TEMPERATURES FOR THREE TESTS 
IN HOT GAS FACILITY. 



(b) MEASURED TRANSIENT LEADING-EDGE THERMOCOUPLE 
DATA. 



(C) SIMULATED THERMAL RESPONSE USED IN STRUCTURAL ANALYSES. 
Figure 2.— Measured and approximated temperature distributions. 
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Figure 3.— Apparent activation energies for creep of copper. Figure 4.— Stress dependence of creep-rate/temperature function c/0. 
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Figure 5.— Elastic analysis of cowl lip, stress in ^-direction. 





Z X 


-.940 -.835 


-.729 -.624 -.519 -.413 

TOTAL STRAIN (X100) 


-.308 


Figure 6 — Elastic analysis of cowl lip, total strain in z-direction. 


-.150 


6 




av 


12.1 


B9 I i 

6.27 .476 


- 5.32 


- 11.1 


STRESS, ksi 


J I I I I L 


- 16.9 


J I L 


- 22.7 


j VO \i 

- 28.5 


83.4 


43.2 


3.3 


- 36.7 - 76.5 

STRESS. MPa 


- 116.5 


- 156.5 


- 196.5 


Figure 7 — Elastic-plastic analysis of cowl lip, stress in z-direction. 
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Figure 8.— Elastic-plastic analysis of cowl lip, total strain in z-direction. 
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Figure 9 — Elastic-plastic-creep analysis of cowl lip, stress in z-direction at 0.75 sec. 
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Figure 10.— Elastic-plastic-creep analysis of cowl lip, total strain in z-direction at 0.75 sec. 
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Figure 1 1.— Elastic-plastic-creep analysis of cowl lip, stress in z-direction at 2.25 sec. 



Figure 12.— Hystersis loops for cowl lip (elastic-plastic-creep analysis, strain 
range Ae = 0.0092). 



Figure 13.— Number of cycles to failure Nf (low-cycle fatigue data obtained 
in argon at 1000 *F (538 # C) with strain rate of 2xl0 -3 sec -1 for 
annealed oxygen-free, high-conductivity copper). 
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